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In this Letter, we study the y-ray signatures subsequent to the production of a Higgs boson in space by dark 
matter annihilations. We investigate the cases where the Higgs boson is produced at rest or slightly boosted and 
show that such configurations can produce characteristic bumps in the y-ray data. These results are relevant in 
the case of the Standard Model-like Higgs boson provided that the dark matter mass is about 63 GeV, 109 GeV 
or 126 GeV, but can be generalised to any other Higgs boson masses. Here, we point out that it may be worth 
looking for a 63 GeV line since it could be the signature of the decay of a Standard Model-like Higgs boson 
produced in space, as in the case of a di-Higgs final state if m% ~ 126 GeV. We show that one can set generic con- 
straints on the Higgs boson production rates using its decay properties. In particular, using the Fermi-L AT data 
from the galactic center, we find that the dark matter annihilation cross section into y+ a Standard Model-like 
Higgs boson produced at rest or near rest cannot exceed (av) ~ a few 10~ 25 cm 3 /s or (av) ~ a few 10~ 27 cm 3 /s 
respectively, providing us with information on the Higgs coupling to the dark matter particle. We conclude that 
Higgs bosons can indeed be used as messengers to explore the dark matter mass range. 



I. INTRODUCTION 

On-going searches at the LHC have been rewarded by one 
of the greatest particle physics discoveries that could possibly 
be made in such a machine, namely the finding of a seem- 
ingly new fundamental scalar or pseudo-scalar particle Orfl]. 
At present measurements of the couplings of this new boson 
to Standard Model (SM) particles along with the absence of 
charged particles tend to suggest that this is a SM Higgs bo- 
son. However this remains to be proven. 

While such a discovery certainly validates our understand- 
ing of the origin of particle masses, it also constrains the types 
of theories that could be proposed to go beyond the Standard 
Model (BSM). For instance, some of the simplest Supersym- 
metric (SUSY) models which have been proposed in the liter- 
ature tend to predict a mass for the Higgs boson that is smaller 
than the measured value m# ~ 125 — 126 GeV |5[] and are 
therefore likely to be ruled out. Moreover, the good agreement 
between the measured branching ratios and those expected in 
the SM (apart perhaps for the two-photon channel) enables 
one to set a stringent constraint on the Higgs invisible decay 
width and to constrain theories in which the Higgs is strongly 
coupled to the dark matter (DM) candidate (%) Ja]. 

Nevertheless, the information collected so far at the LHC 
is not sufficient to exclude the possibility that this new boson 
has a BSM origin. In fact, some non-minimalistic SUSY ex- 
tensions were shown to predict a 'light' Higgs boson with es- 
sentially indistinguishable characteristics from those expected 
within the SM (the remainder of the spectrum in this model 
being typically beyond the scale accessible at LHC) 10]. 
Hence, at present the origin of this new boson remains an open 
question and one needs more clues to determine whether this 
Higgs boson candidate has a SM origin or not. Examining its 



'dark' coupling using other tools than the LHC could be one 
way to proceed. 

In this Letter, we propose to exploit this discovery together 
with recent astrophysical data to constrain the Higgs boson 
production cross section in some specific annihilating DM 
scenarios. We shall focus on the SM-like boson with a mass 
of 126 GeV, but our analysis can be extended to any Higgs 
boson candidate. Now that a SM-like Higgs (or a new) bo- 
son has been discovered and its main characteristics are well 
determined, one can make use of its decay properties (and in 
particular the photon spectrum subsequent to the Higgs bo- 
son decay) to determine whether it has been produced by DM 
in our galactic halo, for instance. Observing the decay of a 
Higgs boson produced at rest (or slightly boosted) in space 
would indeed be suggestive of new physics and provide a new 
window on long-lived neutral particles. The scheme that we 
have in mind is the production of one or two Higgs bosons 
by DM annihilations, although an analogous exercise can be 
done for decaying DM, with similar qualitative arguments for 
DM masses a factor of 2 higher. Once a Higgs boson is pro- 
duced, it is expected to decay immediately, thereby generating 
y-rays. If the associated flux is large enough, this could lead 
to anomalous features in the y-ray spectrum (in particular, an 
excess of photons at some specific energies with respect to the 
background expectations) which can be searched for. Note 
that in what follows we will only focus on the y-ray emission 
from the galactic centre, but our analysis could be extended to 
other regions of the Milky Way as well as the emission arising 
from DM annihilations in dwarf galaxies. 

The y-ray signature associated with a SM-like Higgs bo- 
son decay in our galaxy is expected to be a smooth con- 
tinuum spectrum due to the Higgs decay into SM particles. 
However, if the Higgs boson is produced at rest, one expects 
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FIG. 1 : Examples of Feynman diagrams associated with relevant pro- 
cess discussed in this paper. 



a monochromatic line at Ey ~ 63 GeV due specifically to 
H — > yy, in addition to the continuum. 

The associated signal in an experiment such as Fermi-LAT 
should be a bump around Ey = m# /2 (that is Ey ~ 63 GeV for 
a SM-like Higgs boson) and possibly a broad y-ray excess at 
lower energies, depending on the ratio between the line and 
the continuum. Here we show that it is worth looking for such 
a line in y-ray data, as it could be a mean to probe specific 
annihilating DM scenarios. In particular, in the case of a SM- 
like Higgs boson, one could probe DM masses of about m x ~ 
63 GeV (for %X -> H% m x ~ 109 GeV (for %% -> HZ) or 
m x ~ 126 GeV (for XX -> HH). 

In Section |ll] we discuss the production of the SM-like 
Higgs boson at rest in DM annihilations. After reviewing the 
possible DM annihilation processes which can create one (or 
two) Higgs boson(s) in the final state, we study the detectabil- 
ity of the signature of a Higgs boson decay with the Large 
Area Telescope (Fermi-LAT) on board the Fermi mission and 
discuss the implications for DM scenarios. We also comment 
on the slightly boosted Higgs boson in Section [Til] and con- 
clude in Section [IV] 



H. HIGGS BOSON PRODUCED AT REST BY DM 
ANNIHILATIONS 

In order to produce a Higgs boson in space and at rest, 
the DM mass and spin must have specific values. Quantita- 
tive statements depend on how many Higgs bosons are pro- 
duced in the final state. In the case of DM annihilations 
into two SM-like Higgs bosons, the DM mass must be about 
m% ~ m# — 126 GeV (regardless of its spin). If on the con- 
trary, DM annihilations produce only one SM-like Higgs bo- 
son plus a photon in the final state, the DM mass must be 
about m z ~ mn/2 ~ 63 GeV (assuming that it has a spin- 1/2 
or spin-1) while it should be about 109 GeV if it produces 
a Higgs boson plus a Z boson in the final state (assuming a 
spin-0,1/2 or spin-1). In what follows, we will focus on these 
three specific cases, as they lead to the production of SM-like 
Higgs bosons at rest but, of course, an analogous analysis can 
be done for heavier (presumably BSM) Higgs bosons. We 
now point out some general Higgs boson production mech- 
anisms which could prevail for DM candidates with a mass 
m z ~ 63 GeV, 109 GeV and 126 GeV. Examples of relevant 
Feynman diagrams are given in Fig.[T] 



A. Production mechanisms for m% ~ 126 GeV 

DM candidates with a mass slightly greater than 126 GeV 
can produce two Higgs bosons at rest or near rest in the final 
state either through box diagrams or, if DM is directly coupled 
to the Higgs, through tree-level process (see Fig . |TJ . 

In a SUSY framework for example, two Higgs bosons can 
be produced via box diagrams involving, e.g., charginos and 
W boson or quarks and squarks from the third generation H . 
Disregarding for the moment the possible velocity-squared 
dependence which arises due to the Majorana nature of the 
neutralino, these diagrams are expected to be relatively sup- 
pressed with respect to other annihilation channels which oc- 
cur at tree-level (such as for example neutralino annihilations 
into bb or via a t-channel sbottom or chargino ex- 

change respectively). However they could still be sizable if 
the Higgs boson has large couplings to the particles in the 
box or if there is a large mass degeneracy between the neu- 
tralino and the chargino ( y ± ) for example (if we consider the 
y± _ box diagram la-fllll). Alternatively, the DM could 
also pair annihilate into two Higgs bosons through a pseudo- 
scalar Higgs boson s-channel exchange. If, in particular, the 
mass of the pseudo-scalar is about twice the DM mass, one 
expects a large resonant interaction and potentially a large di- 
Higgs boson production. 

In both cases however, one also expects a large DM pair an- 
nihilation rate into two yy, ZZ, Zy, Hy, HZ leading to extra 
y-ray lines. In many scenarios, these process are related, thus 
giving interesting constraints on the model. However, large 
branching ratios into yy, ZZ, Zy, Hy, HZ could be detrimen- 
tal to the searches for a 63 GeV line. For example, in 'conven- 
tional' BSM scenarios such as SUSY, the di-photon final state 
is supposed to be slightly larger than the di-Higgs produc- 
tion (notably because it is not phase-space suppressed). Since 
the di-photon final state relies on charged loop diagrams, one 
therefore expects a large production of charged particles from 
the DM pair annihilations at tree-level which poses a problem 
for the detectability of the 63 GeV line. Indeed, if the con- 
tribution from annihilations into ^-quarks is significant, it is 
likely that the line at 63 GeV would be totally swamped by 
the continuum y-ray emission resulting from the b hadroniza- 
tion, fragmentation and subsequent decay, with an endpoint 
energy equal to the DM mass, m x ~ 126 GeV. 

There are several ways out, nevertheless. For example, if 
the charged particles which contribute to the direct photon 
emission (loop-suppressed) are all heavier than the DM |[l2tl . 
the DM pair annihilation into such particles is not kinemati- 
cally allowed, thus enabling the di-Higgs final state to be visi- 
ble. In SUSY, this means that one would have to suppress the 
t-channel sbottom exchange diagram and perhaps introduce a 
singlet-like heavy Higgs boson mostly coupled to very heavy 
charged particles 11211 . Alternatively, there could be scenarios 
where the di-photon and di-Higgs final states are produced by 
enhanced box diagrams but in which the sbottom exchanges 
are very suppressed so that the production of fo-quarks at tree- 
level is suppressed. In scenarios with a SM-like Higgs boson 
and no extra pseudo-scalar boson, the tree-level production of 
fo-quarks is expected to be velocity-suppressed. If potential 
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loop/box process, susceptible to imply fe-quarks at tree-level, 
are also suppressed by the introduction of very heavier me- 
diators, the detectability of the 63 GeV line originating from 
enhanced box diagrams could be significant. 

We also note that in models such as the NMSSM where 
one can have both a very heavy (A) and very light (a) pseudo- 
scalar Higgs bosons, the requirement of having a resonant A 
exchange if m z ~ 126 GeV (i.e., niA = 2m% ~ 252 GeV) im- 
plies that one could also produce at tree-level the Aa final 
state, with A produced at rest. The decay of the A into two 
photons could then generate a line at 126 GeV which could be 
confused with the direct (resonant) DM pair annihilations into 
two photons. The dominance of one process over the other 
would mostly depend on the mass difference \mA — 2my\ and 
the strength of the coupling of the neutralino to the Higgs bo- 
son, which itself is constrained by the width of the invisible 
Higgs decay channel lfl3T[l5ll . Such an ambiguity in the ori- 
gin of a possible line at E ~ 126 GeV in this framework could 
be of interest in the context of the 130 GeV and 111 GeV 
bumps observed in the Fermi-LAT data lfl6H2lll . 

For candidates with this mass {m x ~ 126 GeV), the condi- 
tion of predicting a 63 GeV line from a SM-like Higgs bo- 
son produced at rest guaranties a final state with two SM-like 
Higgs bosons. However should such a line be seen, one would 
have to disentangle it from the direct annihilations of DM par- 
ticles with a mass of m x ~ 63 GeV into two photons. Also 
it may be challenging to disentangle the di-Higgs boson final 
state from the Hy final state. These aspects will be discussed 
in the next section. 

Note that all the final states mentioned above have already 
been considered in detail in the literature for generic DM 
masses (see, e.g., Refs. Hi, IH). However, to our knowl- 
edge, the y-ray signature expected from a Higgs boson de- 
cay produced by a ~ 126 GeV DM candidate has not been 
studied explicitly 1 . Many authors have exploited the presence 
of a single photon in DM pair annihilation final states as a 
y-ray signature ifiol I22I I24l - l30ll . However, the possibility of 
these prompt photon lines being accompanied by additional 
lines due to Higgs production at rest has not been pointed out. 
To our knowledge, the fact that the DM pair annihilation into 
two photons could be simply confused with a Higgs boson 
(not necessarily SM-like) production, when ra# ~ 2m x , has 
not been mentioned in the literature yet. 



B. Production mechanisms for m x ~ 63 GeV 

Due to their mass, candidates with m x ~ 63 GeV can only 
produce one SM-like Higgs boson at rest in the final state. 
The DM spin is then fixed by the nature of the second particle 
in the final state. The exact final state can also enable one to 
determine the Higgs boson production mechanism. For exam- 



y-ray fluxes have been calculated and can be obtained from Ref. 12311 for 
different DM masses for DM annihilations into a pair of SM-like Higgs 
bosons, though. 



pie, the Hy final state implies that the Higgs boson production 
must be a loop-suppressed process since the DM is assumed 
to be neutral and cannot produce a photon in the final state 
without coup ling to charged particles (unless one considers 
'dipole' DM (HI). 

Usually one exploits the presence of a single photon in the 
final state to look for such a process (see, e.g., Ref. l22tl ). 
However, the corresponding direct y-ray line would appear at 
very low energy, namely Ey = (1 — m 2 H / (4m?.)) <C 1 GeV, 
to which Fermi-LAT is not sensitive. Hence the only line that 
is experimentally accessible comes from the Higgs decay at 
63 GeV. Nevertheless, observing such a line may not unam- 
biguously point towards the production of a Higgs boson: DM 
pair annihilations into yy could also produce a monochro- 
matic line at the same energy as the Higgs boson decay if the 
DM mass is about 63 GeV. Hence, there could be some confu- 
sion about the origin of the line, even though such a detection 
would definitely point towards new physics. 

In some models, this possible confusion could be solved 
by simply comparing the expected cross sections in different 
channels. For example, in scenarios with photon mixing 1I32I1 . 
the Zd s-channel exchange into yH would be larger than the yy 
final state, so a signal at 63 GeV could be expected. However 
there could be tricky situations. For example, if m x ~ 63 GeV, 
both the XX — > yy and XX Hy process are expected to 
be very large if they are realized through a Higgs portal, i.e., 
XX ~ > H — > yy,Hy. The kinematic condition to see a line at 
m x ~ ma/2 ~ 63 GeV indeed immediately implies that the 
H exchange is resonant. Hence, both final states should be 
copiously produced. If H is the SM Higgs boson, the mag- 
nitude of XX YY versus XX — > Hy is fixed by the ratio of 
the t — t — y versus the t — t —H couplings and the phase space 
factor. Thus, for a SM-like Higgs boson produced very close 
to rest, the phase space factor eventually suppresses a bit the 
Hy final state. Yet, ultimately one should detect the sum of 
the two contributions. 

Note that the importance of the XX ~> YY ar, d XX ^ Hy 
processes through the SM-like Higgs portal ultimately de- 
pends on the mass difference A = 2m x — ma, as well as the 
X — X ~ H coupling. The latter can be tuned (in fact reduced) 
to compensate for the smallness of A, in order to avoid too 
large a resonant annihilation effect, although it cannot be arbi- 
trarily large. The maximum value of the XX^H^yy cross 
section is actually set indirectly by the ATLAS and CMS ex- 
periments. The associated cross section is maximal when both 
A becomes smaller than the Higgs boson decay width (Th) 
and the % — x — H coupling is maximal. Both are being mea- 
sured at LHC through the Higgs visible and invisible decay 
width [33]. A too large X~X~H coupling would make the 
Higgs decay invisible and be in conflict with SM predictions. 

The above discussion assumes that the DM pair annihila- 
tion through the Higgs portal cross section is not velocity- 
dependent. However, if they turn out to be suppressed and 
box diagrams are more important, models with kinetic mixing 
might again lead to a larger value of the cross section for the 
Hy final state (with respect to the yy final state). 

Would such a line be seen, it would remain to be determined 
whether it originates from a SM-like Higgs boson decay into 
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two photons or a model of the type discussed above. However, 
when m x ~ 63 GeV, the DM pair annihilations into any other 
channel would produce a y-ray spectrum with energies Ey < 
m%. Hence the line at ~ 63 GeV would not be buried under 
the continuum spectrum unlike what could occur for m x > 
63 GeV, as discussed in the previous subsection. 



C. Production mechanisms for m % ~ 109 GeV 

When the DM mass is about 109 GeV, the %% -> HZ pro- 
cess can occur (for both bosonic and fermionic DM) via a t- 
channel DM exchange diagram (if DM can couple directly to 
the Higgs) or a s-channel Z exchange diagram. This process 
can also occur through box diagrams. 

For such a value of the DM mass, both the Z and SM-like H 
bosons are produced close to rest and should lead to distinctive 
signatures. In addition to the 63 GeV line from the SM-like 
Higgs boson decay, there could be a line at ~ 109 GeV com- 
ing from the DM annihilations into two photons. Associated 
with this case, there could also be a line at ~ 72 GeV from the 
direct photon in the Hy final state if this channel is not sup- 
pressed. The dominance of one over the other one depends 
again on the couplings and exact process, while their visibil- 
ity essentially depends on the background at these energies. 
Note that y-ray line at ~ 109 GeV from direct annihilation 
into two photons could be consistent with the possible line 
detected at 1 1 1 GeV U3, EH] an d could be used to constrain 
the DM interactions. 



D. Additional remarks 

The results displayed in the next section hold independently 
of whether the new particle discovered at CERN is the Higgs 
boson or not. Since the observed branching ratios are compat- 
ible with the SM Higgs predictions (within 2 a), our conclu- 
sion regarding whether one can see a monochromatic line at 
~ 63 GeV should remain identical. 

Some of the Higgs production mechanisms that we discuss 
in this paper may be associated with a large spin-independent 
elastic scattering cross section with a nucleon and could be 
ruled out by DM direct detection experiments. In particu- 
lar if the DM has a mass in the GeV-TeV range, its interac- 
tions could be severely constrained by the XENON 100 exper- 
iment II341, l35ll . Since this requires to specify a model and we 
intend to set model-independent constraints, we assume that 
the underlying DM particle model is compatible with the re- 
sults from the latest direct detection experiments. However, 
for concrete models such a compatibility has to be checked. 



E. Detectability of the line emission and continuum 

The y-ray emission subsequent to Higgs production typi- 
cally occurs from the Higgs boson decay into, e.g., yy, bb, 
etc. Since all the channels have very well-known branching 
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FIG. 2: Number of photons per GeV produced by the decay of a SM- 
like Higgs boson with mass of 126 GeV produced at rest. The line at 
63 GeV from H — > yy is suppressed, but nevertheless distinguishable 
from the continuum that arises from the Higgs boson decay into the 
rest of SM particles. Note that the very small excess at 30 GeV is 
due to the prompt photon coming from H — > Zy. 



ratios, the y-ray flux can be predicted quite accurately (albeit 
astrophysical uncertainties). 

Predictions depend on the photon energy spectrum 
dNy/dEy associated with the Higgs boson decay. Typically, 
for a Higgs boson of about 126 GeV produced at rest, one 
expects a smooth spectrum (due to dominant decay into b b) 
plus a monochromatic line due to H —> yy IB6T438I1 . In the 
SM, (for niH = 126 GeV) the Higgs boson decay into yy is 
suppressed by a factor of ~ 4 x 10~ 3 with respect to the bb 
final state IB8I1 . so one may think that the y-ray line is hidden 
by the continuum. However, channels such as bb emit pho- 
tons at lower energies than E = mn /2 (owing to final state 
radiation, hadronization, fragmentation and decay). As a re- 
sult, even though the flux associated with the monochromatic 
line is meant to be suppressed, in principle it could be distin- 
guishable from the continuum emission. In order to compute 
the dNy/dEy spectrum, we use PYTHIA 6.4 jH], where we 
set the branching ratio for H -> yy to 2.28 x 10~ 3 Q 2 - The 
result is displayed in Fig. [2] Clearly, the monochromatic line 
appears to be distinguishable from the smooth spectrum, even 
though it is suppressed. 

Now, we estimate the associated flux from DM annihila- 
tions (an analogous analysis could be performed for decay- 



2 Note that the default value in PYTHIA 6.4 is 3.45 x 1CT 3 . Also notice that 
the best fit for the LHC measurement is 2.9 times the SM value and, at 
95% confidence level, it could be up to 5.4 times the SM one |4(J (see also 

Refs. mmi). 
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ing DM) around the galactic center and compare it to the cur- 
rent Fermi-LAT data. We will assume a generic DM candi- 
date, with a thermal average of the annihilation cross section 
times the relative velocity of (av) = (av DMDM ^ H+ ^ z ,m ) = 
3 x 10~ 26 cm 3 /s, where in each case we consider that the only 
annihilation channel is HH, Hy or HZ. 

The differential flux of prompt y-rays generated from DM 
annihilations in the smooth DM halo from a direction within 
a solid angle A£2 is given by lEBTl 

d<t> y (av) dN y If r , , s ,9 

dEy m£ dEy 6K Jaq. J\ s 

where dNy/ dEy is the differential y-ray yield, r) is a symme- 
try factor which for Majorana DM is equal to 1 and 1/2 if 



DM is not a self-conjugate particle, p(r) is the DM density 
profile and r is the distance from the galactic center. The spa- 
tial integration of the square of the DM density profile is per- 
formed along the line of sight within the solid angle of obser- 
vation AO.. More precisely, r = J Rr., — 2sRq cos\|/ + s 2 , and 

the upper limit of integration is s max = J (/^^w — sm2 V^©) + 
i?0cos\|/, where \|/ is the angle between the direction of the 
galactic center and that of observation and Rq is the distance 
from the Sun to the galactic center. Being the contributions 
at large scales negligible, the choice of the size of the Milky 
Way halo, Rmw is not crucial. 

Thus, the flux of DM annihilations can be written as 



d<&y 

dEy 



.27- lO^cm^s -1 XT) x 



dN- 



y {JJ(y)dn\ 



dEy V 20.5 sr / 



100 GeV 

r 



(av) 



Pg 



3- 10- 26 cm 3 /s/ V0.386GeV/cm 3 y V8-25kpc 



with the dimensionless quantity 7(v|/) defined as 



1 



RqP% 



Los 



p( r ( s ,o.)) 2 ds 



(3) 



where for the distance from the Sun to the galactic center and 
for the local DM density we use Rq = 8.25 kpc and p = 
0.386 GeV/cm 3 , respectively JH. 

Although for some DM density profiles, the integration of 
in the solid angle of observation can be done analyti- 
cally j46ll . here we consider an Einasto profile ll47tl . for which 
there is no analytical solution, and compute it numerically. 
This density profile is parametrized as 



p(r) = 0.193 p exp 



- - -1 



oc = 0.17 



(4) 

where r s = 20 kpc is a characteristic length. 

Following Refs. 148145 211 . we consider a 20° x 20° squared 
region centred on the galactic center, for which / J(\\t)dD. = 
20.5 sr. In Fig. [3] we compare the expected flux from this re- 
gion and compare it with the Fermi-LAT data. To obtain the 
measured flux, we take the Fenm'-LAT data obtained from 
August 4, 2008 to October 1, 2012. We extract the data from 
the Fermi Science Support Center archive ll53tl and select only 
events classified as CLEAN. We use a zenith angle cut of 105° 
to avoid contamination by the Earth's albedo and the instru- 
ment response function P7CLEAKLV6. 

In the upper panel of Fig. [3] we show the y-ray spectra for 
three different annihilation channels, Hy (upper red line), HZ 
(black dotted line) and HH (orange line), in which the Higgs 
is produced very close to rest. The DM mass for each case 
is m z = 63 GeV, 109 GeV and 126 GeV, respectively. As 
can be seen from the plot, the fluxes for the three cases are 
very similar, but the Hy final state is slightly more visible 



than the two others 3 , mainly because of the lower value of the 
DM mass in this case. Since the flux scales linearly with the 
cross section, these lines emerge from the y-ray background 
when the associated production cross section is greater than 
(av) ~2.5(5) x 10~ 25 cm 3 /s for Hy (HH), thereby ruling out 
a Higgs boson production cross section larger than this value. 
This can be seen from the lower panel of Fig. [3] where we 
show the value of (a v) for which the signal would be equal to 
the observed background. Interestingly enough, for the case 
of DM annihilations into Hy or HH, producing Higgs at rest, 
the y-ray line from the very suppressed H — > yy channel (see 
Fig. [2j, is expected to provide a more restrictive limit than the 
dominant continuum. Note that this very rough estimate is 
in agreement with the detailed Ferau'-LAT searches of y-ray 
lines 115411 . The limits obtained by the Fermi-LAT analysis for 
m z ~ 63 GeV and %% — > yy can be directly compared to the 
ones presented here for %% — > Hy and m x ~ 63 GeV. While 
the Fermi-LAT limit is (ov) - 3 x 10" 28 cm 3 /s (cf. Fig. 15 
in Ref. JH), we obtain (av) - 2.5 x 10" 25 cm 3 /s, the - 10~ 3 
difference coming from the branching ratio for H — > yy. Sim- 
ilarly, for the case of — > HH and m x ~ 126 GeV, the limit 
obtained from the y-ray line from Higgs decay is just a factor 
of 2 weaker than that for — > Hy and m z ~ 63 GeV (ex- 
plained as a factor of 2 in favour of HH due to having two 
Higgs bosons and a factor of 4 in favour of Hy due to the 
factor of two in the DM mass). 

In DM models where there is a correlation between the 
di-photon and Hy, HZ and/or HH final states, the ratio of 



3 Let us stress again that the y-rays in this case are only those coming from 
Higgs decay. For my ~ 63 GeV, the direct photon lies at energies well 
below detection threshold for Fermi-LAT. 



6 




1 10 

Ey [GeV] 



FIG. 3: Upper panel: Potential y-ray flux from the galactic center 
due to DM annihilating into H y (upper red line), HZ (black dotted 
line) and HH (orange line), when the Higgs is produced very close 
to rest, i.e., for = 63 GeV, 109 GeV and 126 GeV, respectively. 
The results are for an Einasto profile for a 20" x 20" squared region 
around the galactic centre and for (ov) = 3 x 10~ 26 cm 3 /s. The dots 
represent the Fermi-LAT data points in that region for about 4 years. 
Lower panel: Values of the annihilation cross section for each case 
for which the signal flux would be equal to the background flux. Note 
that the astrophysical sources are not subtracted from the data points. 



the flux associated with the prompt y-ray line to that of the 
Higgs boson decay line can be used to test the model. In par- 
ticular when m t ~ 126 GeV, one expects the following ratio 

■p^- = 5=-! — x JUJL Discovering a line at 63 GeV brighter 

than the 126 GeV line would then favour scenarios in which 
the DM pair annihilations into two photons are suppressed by 
at least a factor ~2x 10~ 3 with respect to the di-Higgs pro- 
duction. 

In the absence of evidence for a specific DM model and 
a precise correlation between these two final states, search- 
ing for the Higgs decay line could allow us to obtain a con- 
straint on the DM-Higgs boson interactions. The main dif- 
ficulty associated with these searches consists in removing 
the astrophysical background sources but these searches are 
worthwhile, as they could reveal new physics and point to- 
wards models with multiple scalar and pseudo-scalar Higgs 
bosons with large DM-Higgs couplings, for example. 



III. BOOSTED HIGGS AND MULTIPLE HIGGS BOSONS 
SCENARIOS 

We can now investigate the case of boosted Higgs produc- 
tion and multiple Higgs scenarios. 

A. Boosted Higgs boson 

The Higgs boson decay line considered in the previous sec- 
tion is now replaced by a broad excess which shows up as a 
less prominent feature. For%%^HH, this box-shaped part of 
the spectrum is a particular case of those studied in Ref. f55ll . 
However, in the cases discussed here, this broad excess is ac- 
companied by a smooth spectrum from the Higgs decay into 
all other possible channels plus a possible line due to prompt 
photon emission in the Hy final state. 

These features are illustrated in Fig. 21 where the y-ray 
spectrum due to Higgs decay for a Higgs boson (ra# = 
126 GeV) produced with an energy Eh — 130 GeV is de- 
picted. Over the continuum from the other Higgs decay chan- 
nels, a bump at ~ 60 GeV, corresponding to the Higgs bo- 
son decay into two photons, can still be distinguished. Below 
10 GeV, the continuum is two orders of magnitude (or more) 
brighter than the line, so the limit on the Higgs boson produc- 
tion, for DM masses for which the Higgs boson is boosted, 
is actually obtained from the continuum rather than from the 
broad excess at Ey ~ 60 GeV. This can be seen in Fig. [5] 
which is analogous to Fig. [3] but now for m z = 81 GeV (Hy), 
111 GeV (HZ) and 130 GeV (HH), such that, for all these 
cases, the produced Higgs has an energy close to 130 GeV. 

For the Hy final state, note that there is a y-ray line emitted 
at 32 GeV, in addition to the box-shaped spectrum at Ey ~ 50- 
80 GeV and the continuum from Higgs decays. This line orig- 
inates from the prompt y in the final state and provides the 
most stringent bound on Higgs boson production cross sec- 
tion. Actually, in the case of yy — > Hy, the prompt y-ray is 
always in the energy window accessible by Fermi-LAT if the 
Higgs is not produced very close to rest. Using the Fermi- 
LAT data for this annihilation channel and for m z ~ 8 1 GeV, 
we obtain a limit of about (ov) <4x 10~ 27 cm 3 /s. This is 
comparable to the y-ray line limits obtained by Fermi-LAT 
for y, y, -» yy with m y ~ 32 GeV, that is (ov) <2x 10~ 28 cm 3 /s 
(cf. Fig. 15 in Ref. 115411 ). after correcting the yy — > Hy cross 
section limit by a factor of (1/2) (32/8 1 ) 2 to account for the 
fact that there is only one prompt photon in the Hy final state 
with respect to yy and that the DM mass is different. 

B. Multiple Higgs bosons scenarios 

In minimal SUSY models, in addition to a SM-like Higgs, 
one expects a heavier CP-even Higgs (H2) and a heavier CP- 
odd Higgs (A). If the heavier CP-odd Higgs boson mass is 
about 2m x , annihilations into yy through CP-odd Higgs por- 
tal could be resonant and produce a line at m z . In fact, this 
process has been proposed to explain the bump at 130 GeV in 
the Fermi-LAT data |[56l - l58tl . In these configurations, the Ay 
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FIG. 4: Number of photons per GeV expected from the decay of a 
boosted SM-like Higgs boson with m# = 126 GeV produced with an 
energy Eh = 130 GeV. As one expects, due to the boost, there is no 
line at 63 GeV from H — > yy, but one can nevertheless see a broad 
(box-shaped) emission. 



and H2Y final states might be possible too, leading to the pro- 
duction of a CP-odd Higgs boson on-shell or slightly boosted 
CP-even H2 if niH 2 — niA- These final states should be slightly 
suppressed with respect to the yy final states due to the phase- 
space suppression factor, but would still contribute to the y- 
ray data at Ey = m%. 

In the NMSSM, final states such as A a and Hxa may 
be possible too, with a a second pseudo-scalar Higgs boson 
which can be light and A,H2 two heavy Higgs bosons |[59ll . 
Such final states could lead to the production of Higgs bosons 
produced at rest when 2m x ~ «M,// 2 + m a and could be reso- 
nant when m a <C ttia- The same process could be in fact rel- 
evant for low DM mass scenarios such as those discussed in 
Ref. H. 



IV. CONCLUSIONS 

In this Letter, we have considered the y-ray signatures from 
the decay of a Higgs boson produced in our galactic halo from 
DM annihilations. We have considered, in particular, the case 
where the Higgs boson is SM-like (with a mass of 126 GeV 
and SM branching ratios) and showed that the Higgs boson 
production cross section for annihilating DM particles with 
masses m z ~ 63 GeV, 109 GeV and 126 GeV (Higgs produced 
very close to rest), cannot exceed (av) ~ few x 10~ 25 cm 3 /s. 
The limit is in fact mostly driven by the y-ray line from 
H — > yy. These results can be trivially generalised to other 
Higgs boson masses (as relevant in BSM models with multiple 
Higgs bosons and Higgs mass spectrum such as the NMSSM) 




1 10 100 

Ey [GeV] 



FIG. 5: Same as Fig.[5] but for DM masses such that the Higgs boson 
is produced with an energy Eh — 130 GeV, i.e., for m x = 81 GeV 
(Hy, upper red line), 1 1 1 GeV (HZ, black dotted line) and 130 GeV 
(HH, orange line). 



leading to different DM scenarios. 

We have also considered the case of a slightly boosted 
Higgs boson and shown that the associated signature would 
exhibit a broad (box-shaped) y-ray excess. However, the con- 
tinuum associated with the other Higgs boson decay modes 
and to the second particle in the final state would lead to a 
brighter y-ray emission, which can be used to constrain the 
Higgs boson production cross section. Focusing in particu- 
lar on the Hy final state for a SM-like Higgs boson produced 
with an energy Eh = 130 GeV, we find that the Higgs boson 
production cross section cannot exceed ~4x 10~ 27 cm 3 /s. 

Therefore, we have obtained a simple estimate for the limit 
on the Higgs boson production cross section that is indepen- 
dent of any other DM annihilation channels and demonstrates 
that performing Higgs boson decay line searches could be use- 
ful to probe the Higgs boson dark couplings (i.e., couplings to 
DM particles). This must be compared to the limits set on the 
invisible Higgs boson decay branching ratios obtained by us- 
ing LHC measurements (cf., for example, Ref. lfl5tl ). but the 
two approaches (collider and indirect detection searches) are 
complementary. 
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